Investigation on nonlinear lateral-torsional coupled vibration of a rotor system with substantial unbalance

Substantial unbalance may be caused by fan blade off during the operation period of gas turbine engines, and related dynamic problems are very critical to the safety design of rotor system in aero-engine. This article aims to understand lateral-torsional coupled vibration of the rotor system with substantial unbalance. The governing equation of a modified unbalanced rotor system is established based on Lagrangian approach. Then, a mathematical analytical method is proposed in which a linear approximation is derived and the Floquet theory and Hill’s method are incorporated, from which the modal characteristics of the unbalanced rotor are obtained. The modal characteristics of the unbalanced rotor system are revealed comprehensively for the first time. Furthermore, the relation between the modes and responses of the unbalanced rotor is discussed in detail. The results show that the lateral vibration and torsional vibration of the unbalanced rotor are coupled through the inertial terms in the governing equations. Due to the coupling, veering and lock-in phenomena occur between the frequencies of the forward whirl mode and the torsional mode. Furthermore, lock-in can lead to a kind of principal instability. With regard to the response of the unbalanced rotor, both natural vibration components and enforced vibration components appear in the lateral response, while only natural vibration components appear during torsional vibration. Moreover, natural vibration components play a crucial role in the response within the principal instability region and cause divergence of the vibration amplitudes in the lateral and torsional directions.     

针栓式推力室冷却特性试验研究

针栓式发动机具有推力调节简单、声学燃烧稳定性好和成本低廉等一系列优点,但目前对针栓式推力室的冷却特性掌握还不够深入,研制过程中曾多次出现因推力室冷却组织不好而烧蚀的问题。采用针栓式推力室试验件,对身部的冷却特性进行了深入研究,结果表明:推力室圆柱段前端的温度较低,无需采用专门的冷却措施;末端内侧气壁温度会达到约1650℃,局部存在明显的烧蚀现象,必须采取有效的热防护措施;下游冷却液膜量的变化对圆柱段的冷却特性影响较小。     

Applications of structural efficiency assessment method on structural-mechanical characteristics integrated design in aero-engines

In the design process of advanced aero-engines, it is necessary to carry out an effective analysis method between structural features and mechanical characteristics for a better structural optimization. Based on the structural composition and functions of aero-engines, the concept and contents of structural efficiency can reflect the relation between structural features and mechanical characteristics. In order to achieve the integrated design of structural and mechanical characteristics, one quantitative analysis method called Structural Efficiency Assessment Method (SEAM) was put forward. The structural efficiency coefficient was obtained by synthesizing the parameters to quantitatively evaluate the aero-engine structure design level. Parameterization method to evaluate structural design quality was realized. After analyzing the structural features of an actual dual-rotor system in typical high bypass ratio turbofan engines, the mechanical characteristics and structural efficiency coefficient were calculated. Structural efficiency coefficient of high-pressure rotor (0.43) is higher than that of low-pressure rotor (0.29), which directly shows the performance of the former is better, there is room for improvement in structural design of the low-pressure rotor. Thus the direction of structural optimization was pointed out. The applications of SEAM shows that the method is operational and effective in the evaluation and improvement of structural design.     

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符；反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。      

弱引力下小行星动量驱动机器人的可控跳跃行为

以小行星表面着陆探测为背景,提出一种动量驱动机器人(MoRo)以满足弱引力复杂环境下的探测需求。该机器人利用弱引力环境下的摩擦和碰撞特性,通过主动辨识环境参数,规划和控制动量轮以产生期望的驱动力矩,完成可控性跳跃及腾空后的稳定拍照等任务。首先,基于MoRo的动量轮刹车机构特性,分析了MoRo在弱引力环境下的跳跃机理并对其跳跃方式进行了规划;接着考虑动量轮驱动机构三闭环伺服系统的非线性特性,基于Herze碰撞模型和Karnopp摩擦模型建立了MoRo在小行星表面的跳跃行为动力学模型;其次,使用机器学习算法建立环境参数和MoRo运动的函数关系,并基于环境参数规划动量轮转速实现跳跃距离和腾空高度的可控。最后,通过数值仿真校验了MoRo跳跃规划方法和控制方法的可行性。     

空天飞行器返回滑翔段在线制导方法

针对空天往返飞行器的返回滑翔段在线制导问题，设计了一种新的滑翔段飞行剖面，实现了滑翔段终端交班高度、位置和倾角约束的自动满足，减少了在线制导算法中需处理的约束数量。推导了滑翔段运动状态、过程约束和性能指标的解析表达式，获得了剩余航程和终端速度间的函数关系。在此基础上，提出了一种双层在线制导方法：内层解析重构飞行剖面，同时通过解析确定路径点来改变剩余航程的变化率，进而对终端交班速度进行控制；外层借助解析表达式，使用粒子群优化算法(PSO)和改进共轭梯度法(CGM)优化飞行剖面，从而满足过程约束和指标要求。最后通过数学仿真验证了方法的正确性。     

全模颤振风洞试验三索悬挂系统多体动力学分析

全模颤振风洞试验需要通过软支撑系统模拟飞行器的自由飞行状态并调整模型姿态达到配平状态。参考NASA双索悬挂方案，提出了一种两电机驱动的三索悬挂系统，利用后方两索的同向/反向联动实现模型俯仰和滚转姿态的调整，利用弹簧刚度以及钢绳张力设计实现支撑频率要求。基于柔性多体动力学方法，建立了包括飞行器刚体模型、柔性索、滑轮、弹簧、气动力模型、伺服电机控制在内的复杂系统动力学模型，其中，利用任意拉格朗日-欧拉（ALE）变长度索单元描述钢绳，利用不约束物质坐标的索结点约束描述钢绳与滑轮相互作用，利用索结点物质输运速度约束描述伺服电机绞盘，利用飞行力学的气动力模型描述吹风下的气动力。基于该模型，通过小扰动响应辨识研究了弹簧刚度、钢绳张力、连接点位置等因素对支撑频率的影响规律，并分析了系统姿态调整能力，俯仰调整范围达到-12.5°~12.5°，滚转调整范围达到-45°~45°。采用滑轮处电位计测量的钢绳相对位移作为反馈信号，基于设计的控制律利用多体动力学求解器与Simulink对风洞吹风下的姿态调整过程进行仿真，模型达到配平状态，获得了吹风下的索拉力和伺服电机功率，为系统设计提供基础。     

基于合成射流的旋翼翼型动态失速控制研究

针对直升机旋翼工作环境下来流速度和迎角(Angle of attack,AoA)耦合引起的动态失速问题,建立了基于合成射流的旋翼动态失速控制的数值模拟方法。采用运动嵌套网格方法,通过对翼型的平移和旋转实现变来流速度-变迎角的耦合。以积分形式的雷诺平均N-S方程为主控方程,空间离散使用Roe格式,时间离散为隐式LU-SGS方法,以OA209翼型为研究对象,在翼型上表面放置合成射流激振器,开展了射流位置、动量系数、无量纲频率以及偏角等参数对轻度失速、深度失速下翼型动态失速控制的研究。研究发现,轻度失速下,射流位置靠近气流分离点时(20%c附近,c为翼型弦长),对逆压梯度引起的轻度失速控制效果最佳。深度失速下气流分离点虽在5%c之前,但射流位于前缘分离泡后端(10%c附近)时控制效果较好。大迎角需要较大的动量系数才能有效控制。射流频率对涡结构的尺寸和数量会产生一定影响,能改变气动特性波动幅度。较小的射流偏角对轻度失速的控制更有效,而深度失速则需要较大的偏角。     

多自由度振动环境下惯组非线性传递特性研究

惯组作为飞行器姿控系统的传感器,其局部安装结构的传递特性的测量精度直接关系到导航精度。目前,惯组普遍使用减振器进行隔振,而减振器都呈现出较强的非线性特征。为了考察惯组在不同工况下的传递特性,将惯组简化为六自由度Duffing模型,推导了基础激励下系统的运动微分方程,并用龙格-库塔法对方程进行求解,分析了自由衰减振动和强迫振动下不同工况的系统传递特性。结果表明,多自由度激励比单自由度激励工况得到的系统传递特性的频率和幅值都低。考虑到惯组真实的使用环境,应当在多自由度振动环境下进行传递特性试验。     

Analysis of the dynamic behavior of structures using the high-rate GNSS-PPP method combined with a wavelet-neural model: Numerical simulation and experimental tests

Recently, the high rate global navigation satellite system-precise point positioning (GNSS-PPP) technique has been used to detect the dynamic behavior of structures. This study aimed to increase the accuracy of the extraction oscillation properties of structural movements based on the high-rate (10?Hz) GNSS-PPP monitoring technique. A developmental model based on the combination of wavelet package transformation (WPT) de-noising and neural network prediction (NN) was proposed to improve the dynamic behavior of structures for GNSS-PPP method. A complicated numerical simulation involving highly noisy data and 13 experimental cases with different loads were utilized to confirm the efficiency of the proposed model design and the monitoring technique in detecting the dynamic behavior of structures. The results revealed that, when combined with the proposed model, GNSS-PPP method can be used to accurately detect the dynamic behavior of engineering structures as an alternative to relative GNSS method.